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Pre-mRNA splicing is required for the accurate expression of virtually all human protein coding genes. However, splicing 
also plays important roles in coordinating subsequent steps of pre-mRNA processing such as polyadenylation and mRNA 
export. Here, we test the hypothesis that nuclear pre-mRNA processing influences the polyribosome association of al- 
ternative mRNA isoforms. By comparing isoform ratios in cytoplasmic and polyribosomal extracts, we determined that 
the alternative products of -30% [597/1954) of mRNA processing events are differentially partitioned between these 
subcellular fractions. Many of the events exhibiting isoform-specific polyribosome association are highly conserved across 
mammalian genomes, underscoring their possible biological importance. We find that differences in polyribosome as- 
sociation may be explained, at least in part by the observation that alternative splicing alters the c/5-reguIatory landscape 
of mRNAs isoforms. For example, inclusion or exclusion of upstream open reading frames (uORFs) in the 5'UTR as well as 
A/w-elements and microRNA target sites in the 3'UTR have a strong influence on polyribosome association of alternative 
mRNA isoforms. Taken together, our data demonstrate for the first time the potential link between alternative splicing 
and translational control of the resultant mRNA isoforms. 



[Supplemental material is available for this article.] 

The coding sequences (exons) of most metazoan genes are sepa- 
rated by intervening sequences (introns) that must be removed 
from the precursor message RNA (pre-mRNA) prior to translation. 
Pre-mRNA splicing is an essential step in gene expression and is 
catalyzed by the spliceosome, a dynamic macromolecular ribo- 
nucleoprotein particle, which assembles de novo on each intron 
(Hoskins and Moore 2012). Pre-mRNA splicing is also a major 
regulatory point in the gene expression pathway (Braunschweig 
et al. 2013). The vast majority of multiexon genes are alternatively 
spliced, giving rise to a diverse but functionally enigmatic pool of 
transcripts (Castle et al. 2008; Pan et al. 2008; Wang et al. 2008). 

In eukaryotic systems, there is a poor correlation between 
steady-state mRNA and protein levels (Lian et al. 2001; Griffin et al. 
2002; Cox et al. 2005; Kislinger et al. 2006; Wilhelm et al. 2006; 
Schmidt et al. 2007; Fu et al. 2009). Although some of this dis- 
crepancy can be attributed to mRNA stability and the efficiency of 
translation initiation, the contribution of alternative splicing (AS) 
to this problem is largely ignored. Conversely, the impact of AS on 
mRNA diversity is indisputable and yet the underlying biological 
relevance for many transcript isoforms remains enigmatic. Bridg- 
ing this gap between transcript diversity and protein complexity 
requires a closer examination of isoform-specific mRNA fates. 

AS regulates the expression of human genes by at least two 
well-established mechanisms. First, AS expands protein diversity 
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by generating different mRNA isoforms from protein-coding genes. 
Second, AS attenuates gene expression by producing substrates 
for nonsense-mediated decay (AS-NMD) (Ni et al. 2007; Pan et al. 
2008; Wang et al. 2008). In contrast to NMD, little is known re- 
garding how mRNA isoforms differ in their translatability. It is clear 
that the trans-acting factors associated with alternative mRNA 
isoforms can vary just as distinctly as the protein-coding capacities 
of the messages they encode. Many RNA binding proteins, in- 
cluding shuttling proteins, are implicated in multiple steps of post- 
transcriptional gene regulation, including mRNA export, stability, 
surveillance, localization, and translation (Sanford and Caceres 
2004; Zhang and Krainer 2004; Zhong et al. 2009; Han et al. 2011). 
These factors can be transcript-associated throughout the life cycle 
of an mRNA and may coordinate nuclear and cytoplasmic steps of 
gene expression (Dreyfuss et al. 2002; Nashchekin et al. 2007; 
Kylberg et al. 2008). This hypothesis is supported by the observa- 
tion that exonic splicing enhancers can also alter the efficiency of 
mRNA translation in a position-specific manner (Sanford et al. 
2004; Michlewski et al. 2008; Delestienne et al. 2010). Thus, an 
unanticipated consequence of alternative splicing may lie within 
the subsequent regulation of the RNA binding proteins packaged 
upon each isoform: functionally distinct mRNPs with distinct cy- 
toplasmic fates (Zhong et al. 2009; Wang et al. 2012). 

Here, we test the hypothesis that alternative mRNA isoforms 
are subject to translational regulation. This approach combines 
subcellular frac tionation and high throughput RNA sequencing 
(Frac-seq) to assay partitioning of mRNA isoforms between the 
cytosolic and the polyribosome-associated fractions. Our data 
shows that —30% of alternative splicing events exhibit isoform- 
specific polyribosome association. The sequences encompassing 
these events are more evolutionarily conserved compared to ar- 
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chetypal alternative pre-mRNA processing events, underscoring 
their biological importance. We also suggest that alternative splicing 
alters the ds-regulatory landscape of mRNA isoforms, leading to 
distinct cytoplasmic fates. Taken together, our results document an 
unexpected link between alternative pre-mRNA processing and 
the cytoplasmic fate of endogenous mRNA isoforms. 

Results 

Frac-seq reveals isoform-specific differences in polyribosome 
association 

Polyribosome profiling is widely used to estimate the translation 
status of mRNAs (Masek et al. 2011). We hypothesized that alter- 
native splicing not only expands the protein coding capacity of 
human genes, but also establishes patterns of isoform-specific 
translational control. To test this hypothesis, we propose a geno- 
mic approach comparing alternative mRNA isoform diversity in 
cytoplasmic and polyribosome subcellular fractions (Fig. 1A). The 



former sample represents the output of the nucleus, whereas the 
latter represents those mRNAs that are engaged in the translation 
machinery. We extracted only the polyribosomal peaks from each 
gradient thereby avoiding possible contaminating RNA from other 
regions of the gradient (Supplemental Fig. 1). 

Following paired-end sequencing of polyA + selected RNA li- 
braries and mapping to hgl9 (Supplemental Table 1), we employed 
a previously developed Bayesian inference methodology to esti- 
mate percentage spliced-in MISO) (Katz et al. 2010) in the 
polyribosomal and cytoplasmic fractions. The value of *P was es- 
timated for eight different types of alternative events, including 
skipped exons (SE), alternative 5' splice sites (A5SS), alternative 
3 'splice sites (A3SS), mutually exclusive exons (MXE), retained 
introns (RI), alternative first exons (AFEs), alternative last exons 
(ALEs), and tandem 3'UTRs (alternative polyadenylation). We 
found that the values estimated for biological replicates are 
strongly correlated, suggesting a sufficient degree of biological re- 
producibility in the data (Supplemental Fig. 2). 
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Figure 1. Frac-seq reveals isoform-specific differences in polyribosome association. (A) Schematic diagram of the Frac-seq protocol. (B) Bayes factor 
threshold estimation as a function of biological reproducibility (measured via Spearman rank correlation). Local smoothing of the mean and standard error 
of the mean (SEM) provided as line and shading, respectively. (C) Empirical false discovery estimation (FDR) using bootstrapped (1 OOx) label permutation 
of one biological replicate for each |A^| cutoff. Mean FDR and standard error are plotted as solid and dashed lines, respectively. (D) Percentage of 
alternative events showing significant changes in polyribosome association (blue) compared to detectable alternatively spliced events showing no dif- 
ference (gray), (f ) Complexity of RNA processing events per gene (columns scaled to percentage), with complex events (more than one differential event 
per gene) in gray, and simple (one differential event per gene) events in blue. (F) Breakdown of significant event types and their distributions with respect 
to A^P. Schematics on the left side present the "spliced in" region as darkly shaded. For example, the skipped exon (SE) diagram presents the shaded exon 
that when preferentially associated with polyribosomes would show a positive value for A^. 



1 61 6 Genome Research 



www.genome.org 



Coupling of alternative splicing and translation 



To determine which events, if any, showed significant changes 
in ^ upon comparison of polyribosomal and cytoplasmic frac- 
tions, we utilized the Bayes factor (BF) as a measure of statistical 
significance for A^ (^ Po iyribosomai - ^cytoplasm) (Fig- IB). As a pre- 
requisite for biological significance, we assume that the events 
must be biologically reproducible. Proceeding from this logic, we 
split the data (BF > 1 in both replicates) into vigintiles following the 
distribution of Bayes factor values. Plotting Spearman's rank cor- 
relation coefficient (A^ replicate I vs. replicate II) for the events in 
each quantile produces a hyperbolic curve that approaches perfect 
correlation and plateaus near a Bayes factor of five (Fig. IB). These 
data suggest that the Bayes factor is indeed a measure that corre- 
sponds to increasing biological reproducibility, and a cutoff of five 
is an appropriate threshold. To determine the rate of false discovery 
(FDR), we utilized the real data as a set of true positives and esti- 
mated false positives from 100X bootstrapped label permutations. 
Introducing |A^| consistency between replicates as a second strin- 
gency measure, we plotted this empirical FDR for each boot- 
strapped replicate following increasing cutoffs for A^P (Fig. 1C; see 
Methods). These data allowed us to choose relevant criteria for 
significance (BF > 5, |A^| > 0.1) that maximize the number of sig- 
nificant events and minimize the false discovery rate to an ac- 
ceptable level (FDR « 5.02 ± 0.9%). 

If mRNA translation is independent of the process of alter- 
native splicing, one prediction would be that the relative fractional 
abundance of alternatively spliced messages would be consistent 
in both the experimental fractions. Using Frac-seq, we identified 
597 alternative events with significant difference in ^ between the 
cytoplasmic and polyribosomal fractions (Supplemental Table 2). 
We also identified 1357 alternative processing events (0.1 < ^ < 
0.9) that are unchanged between the two fractions ("no-diff", |A^| 
< 0.1; BF < 1). Taken together, we find that 30.5% of detectable 
alternative processing events are differentially partitioned between 



the cytosolic and polyribosomal fractions (Fig. ID). Further in- 
vestigation showed that >40% of these events are the only differ- 
ential event per gene (Fig. IE, "simple events"). This observation 
suggests that loss or gain of sequence elements associated with 
single alternative events may contribute to the isoform-specific 
polyribosome association we detect. The 597 significant events 
encompass all eight distinct pre-mRNA processing classes defined 
in the MISO database. Furthermore, we recognize a marked pref- 
erence for significant changes in ^ affiliated with alternative last 
exon (ALE) and alternative first exon (AFE) events (Fig. IF; Sup- 
plemental Fig. 3). This is consistent with roles for 5' and 3' UTRs 
(untranslated regions) in modulation of both mRNA stability and 
translation. Additionally, tandem UTR events (alternative poly- 
adenylation) showed a clear bias toward increased polyribosome 
association for the shorter isoform. This observation is consistent 
with previous reports (Sandberg et al. 2008; Mayr and Bartel 2009). 
Collectively, these data strongly support our hypothesis that al- 
ternative pre-mRNA splicing impacts the translation of mRNA. 
Finally, genes associated with significant differences in isoform- 
specific polyrbosome association are enriched in gene ontology 
terms related to RNA binding and pre-mRNA splicing (Supple- 
mental Table 4). 

Validation of alternative event-specific polyribosome 
association by RT-PCR 

We tested seven alternative cassette exons and one alternative last 
exon (ALE) estimated by MISO to be differentially partitioned be- 
tween the cytoplasmic and polyribosomal fractions (Supplemental 
Fig. 4). Qualitatively, the RT-PCR analysis for seven of the eight 
events was in agreement with the computational estimates (sum- 
marized in Supplemental Table 3; Supplemental Fig. 4). Figure 2A 
compares A^ estimates from sequencing reads and RT-PCR for 
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Figure 2. Validation of alternative mRNA isoforms in polyribosomal and cytosolic subcellular fractions by RT-PCR. (A) Correlation of endpoint RT-PCR (in 
triplicate with standard error) determined A^ with average MISO estimated A^. Linear regression plotted as a solid line with 95% confidence interval 
(shaded region). Significance P-value derived from t-distribution test using Pearson's product-moment correlation coefficient. (B) Representative agarose 
gels of endpoint RT-PCR for two exemplar genes, HNRNPR and KARS. CI and C2 refer to the cytoplasmic derived replicates, whereas PI and P2 refer to the 
polyribosome-bound replicates. (C,D) Sashimi plots depict RPKM coverage of RNA-seq reads used to estimate M> values for HNRNPR and KARS, respectively. 
Junction reads are annotated as connecting lines. 
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seven cassette exons. This provides a —12.5% experimental esti- 
mate of the FDR that is not unexpected given the sample size and 
our computational estimate. One splicing event in the MOV10 
gene was clearly a false positive and another, in the SGOL1 gene, 
proved to be qualitatively similar to the MISO estimate but was 
difficult to quantify by RT-PCR (Supplemental Fig. 4). Despite these 
two caveats and a relatively small sample size, we observe signifi- 
cant positive correlation between RT-PCR and computational es- 
timates of (Fig. 2 A; Supplemental Table 3). Two exemplars, 
HNRNPR and KARS, are shown in Figure 2B and represent robust 
examples of event-specific polyribosome association. Computa- 
tional estimates of ^ suggest that 15%-32% of HNRNPR transcripts 
contain the alternative exon in the polyribosome fraction com- 
pared to 58%-66% in the cytoplasm (Fig. 2C). Alternative mRNA 
isoforms produced by the KARS gene show the opposite pattern of 
partitioning. Estimates of *P suggest that 46%-64% of KARS poly- 
ribosome associated transcripts contain the alternative cassette 
exon compared to only 23%-31% in the cytoplasm (Fig. 2D). 
Representative endpoint RT-PCR confirmed these results (Fig. 2B). 

Translational control can be assayed by determining the dis- 
tribution of mRNAs across individual polyribosome fractions. One 
caveat of Frac-seq is that the polyribosome fractions are pooled 
prior to sequencing. To determine if isoform-specific polyribosome 
association is indicative of translational control, we examined the 
sucrose gradient distributions of alternative mRNAs from four 
genes: DDX50, DGUOK, HNRNPR, and KARS. Quantification of RT- 
PCR analysis of skipping events in different fractions demonstrates 
that polyribosome association of alternative mRNA isoforms is 
highly variable (Supplemental Fig. 5). For example, the DGUOK 
isoform including a cassette exon is associated almost exclusively 
with the heavier polyribosome fractions (Fig. 3, lanes 5-8). In 
contrast, the isoform lacking this exon is distributed fairly evenly 
across the gradient, suggesting that the skipped isoform is trans- 
lated less efficiently than the included isoform. Contrary to DGUOK, 
both KARS and DDX50 behave similarly with increased exon in- 
clusion in the lighter polyribosome fractions and enrichment of 
the skipped isoform in the heavier fractions (Fig. 3, lanes 1-4 and 
5-8, respectively). For HNRNPR, inclusion of the alternative cas- 
sette exon is enriched in the cytoplasmic fraction and is dramati- 
cally reduced in the polyribosome fractions (Fig. 3, cf. lane 9 to 
lanes 1-7). 

Sequences corresponding to isoform-specific polyribosome 
association are conserved through evolution 

Although our data demonstrate that —30% of alternatively pro- 
cessed messages are differentially associated with polyribosomes, 
the mechanisms by which they do so remain elusive. Phylogenic 
conservation of DNA sequences often underscores their functional 
importance. This hallmark has been successfully applied to alter- 
native splicing to identify putative regulatory elements in alter- 
native exons and flanking intronic sequences (Sugnet et al. 2006; 
Wang et al. 2008). We compared the PhyloP conservation scores 
(placental mammals using a 46-way alignment) for different clas- 
ses of alternative events associated with isoform-specific poly- 
ribosome association (Siepel et al. 2005; Pollard et al. 2010). As 
expected, cassette exons appear to be less conserved than adjacent 
constitutive exons, suggesting recent evolution (Gelfman et al. 
2012). Surprisingly, we observe elevated conservation scores for 
differential cassette exons in the coding sequences flanking the 
upstream intron (Fig. 4A, gold and blue lines, respectively). Most 
notably, we detect convincing changes in conservation near the 
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Figure 3. Distribution of alternative mRNA isoforms in ranked poly- 
ribosome complexes. (Upper panel) Absorbance profile (A 2 6o) of cytosolic 
extracts across 15%-45% sucrose gradient. Purified polyribosome frac- 
tions are indicated by P2 to >P9 corresponding to two ribosomes (P2) up 
to nine or more (^P9). (Lower panels) Alternative splicing of KARS, 
HNRNPR, DGUOK, and DDX50 in different polyribosomal peak fractions 
assayed by RT-PCR and agarose gel electrophoresis. The constitutively 
spliced ACTB transcript is present at the bottom. 

transcription start site (TSS) for the proximal alternative first exon 
and 250-bp internal to the upstream alternative polyadenylation 
(PA) site (Fig. 4B,C). In both cases, sequences enriched in the 
polyribosome fraction have higher conservation scores compared 
to those in the "no-diff ' set (Fig. 4B,C, respectively). In contrast, we 
observe sub-background levels of conservation present in the 
proximal ALE and subtle increases in conservation 300 bp up- 
stream of the PA site in the distal ALE (Fig. 4D). Taken together, our 
data suggest that splicing events related to strong polyribosome 
association show evolutionary conservation, which may reflect 
the biological importance and/or functionality of those sequences. 

Alternative pre-mRNA processing alters the c/5-reguIatory 
landscape of mRNA isoforms 

A significant fraction of alternative splicing events give rise to 
nonfunctional mRNAs containing premature termination codons 
(PTC). These potentially deleterious mRNA isoforms are highly 
unstable due to their efficient elimination by the nonsense- 
mediated mRNA decay pathway. Transcripts containing PTCs are 
hypothesized to be degraded in a pioneer round of mRNA trans- 
lation and thus should be largely excluded from the polyribosome 
fraction. To identify possible NMD substrates, we searched for al- 
ternative splicing events that introduced sequences containing 
termination codons in all three frames. We observe that alternative 
events containing PTCs are depleted from the polyribosome frac- 
tion (Fig. 5). Given their prominent role in translational control 
and mRNA stability, we also analyzed alternative events that remodel 
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Figure 4. Conservation of DNA sequences encompassing different classes of polyribosome associ- 
ated alternative mRNA processing events. Local nonlinear smoothing using a generalized additive 
model shows best-fit sequence conservation (46-way PhyloP) scores of alternative events and flanking 
sequence. Significant events are depicted as yellow, while "no diff" events are colored blue (shaded 
95% confidence interval). (A) Skipped exon events; (B) alternative first exon events; (C) alternative last 
exon events; and (D) alternative polyadenylation events. 



sequences in 5' and 3' UTRs associated with isoform-specific 
polyribosome association. We exploited a list of predicted func- 
tional elements annotated in the UTR database (UTRdb) (Grillo 
et al. 2009) and intersected these annotations with the set of al- 
ternative events exhibiting a change in polyribosome association 
(see Methods). Included in this list are well-annotated 5' and 3' 



UTR specific regulatory elements such 
as upstream ORFs (uORFs) and internal 
ribosomal entry site (IRES), as well as 
microRNA target sites and Alu elements, 
respectively. These element classes are 
well characterized for their roles in trans- 
lational control and mRNA stability (Grillo 
et al. 2009; Gong and Maquat 2011a,b; 
Ramani et al. 2011). Our analysis dem- 
onstrates that uORFs and IRES elements 
are associated with an increase in polyri- 
bosome association for element-harboring 
isoforms compared to the element-de- 
ficient isoform (P- value < 3.79 X 10~ 3 and 
<1.15 X 10~ 3 , respectively) (Fig. 5). For 
alternative pre-mRNA processing events 
within the 3 'UTR, we observe a signifi- 
cant decrease in polyribosome associa- 
tion for isoforms containing microRNA 
target sites or Alu elements when com- 
pared to the isoform lacking these ele- 
ments (P-value < 6.42 X 10" 6 and <6.34 X 
10~ 3 , respectively) (Fig. 5). 

The analysis described above suggests 
that remodeling of known ds-regulatory 
elements may contribute to event-specific 
patterns of polyribosome association. To 
identify novel ds-elements, we searched 
sequences associated with preferential 
polyribosome or cytoplasmic distribution 
for over-represented n-mers. It was most 
feasible to analyze the ALE and AFE events 
given the relatively small sample size of 
other event types. We identified several 
distinct hexamer sequence clusters that 
are over-represented in either polyri- 
bosome or cytoplasmic fractions for AFE 
events (Supplemental Fig. 6A). These se- 
quences correspond primarily to GC-rich 
hexamers. The only enriched hexamers we 
identified in cytoplasmic or polyribosomal- 
associated ALE events resembled the 
ds-elements related to polyadenylation 
signals (Supplemental Fig. 6C). 

Discussion 

Here we present Frac-seq, a method to 
determine the enrichment or depletion 
of alternatively spliced mRNAs in poly- 
ribosomes. The Frac-seq method is related 
to the recently reported ribosome foot- 
printing method but differs in several im- 
portant ways (Ingolia et al. 2009; Ingolia 
2010). First, Frac-seq is not intended to 
quantify the translational state of all 
mRNAs but rather to determine the alternatively spliced tran- 
scriptome that is associated with the translation machinery. The 
footprinting method produces a ribosomal silhouette within the 
coding sequence of each mRNA that provides an extremely accu- 
rate measure of ribosome density. However the nuclease digestion 
decreases the complexity of the library and lowers the resolution of 
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Figure 5. Alternative pre-mRNA processing alters the c/s-regulatory 
landscape of mRNA isoforms. Intersection of the alternative portion of each 
significant event with predicted c/s-acting RNA elements. Events whose al- 
ternative region harbors or overlaps a c/s-acting element are labeled by their 
direction of change (positive or negative A^P, polyribosomal or cytoplasmic, 
respectively). Premature termination codon (PTC) containing events have at 
least one termination codon in all three reading frames of the alternative 
portion of the event. Events with uORF annotations contain the start codon of 
the annotated upstream open reading frame. The tendency of an annotated 
element-containing sequence toward the cytoplasm is colored gray, while 
the polyribosome is dark gray. Results of a x goodness-of-fit test to an even 
ratio are denoted as follows: (***) <0.001; (**) <0.01; (*) <0.05. 

different mRNA isoforms. Additionally, information pertaining to 
processing within UTRs is missing, as these regions are not pro- 
tected from nuclease. 

In eukaryotic systems, mRNA transcript levels often do not 
correspond with protein quantities (Lian et al. 2001; Griffin et al. 
2002; Cox et al. 2005; Kislinger et al. 2006; Wilhelm et al. 2006; 
Schmidt et al. 2007; Fu et al. 2009). For example, in human me- 
dulloblastoma cells, there is a relatively poor correlation between 
steady-state protein and mRNA levels (R = 0.29) (Vogel et al. 2010). 
The authors of this study identified sequence features such as 
transcript length, secondary structure and known regulatory ele- 
ments accounting for up to 67% of the variation between mRNA 
and protein levels (Vogel et al. 2010). However the contribution of 
alternative splicing (AS) to this problem is poorly understood. 
Vogel et al. (2010) demonstrate that coding and 3'UTR sequences 
explain a significant fraction of protein abundance variation. 
Our data largely agree with these results. In particular, we find that 
isoforms enriched in the polyribosome fraction frequently corre- 
spond to exon skipping events or processing events leading to 
distinct or shorter 3 'UTRs (Fig. IF). However, because we adopted 
an event-based analysis, our results do not address the relationship 
between isoform length and polyribosome association. Based on 
results from ribosome footprinting experiments in yeast, we spec- 
ulate that isoforms with shorter coding regions will be translated 
more efficiently (Ingolia et al. 2009). 

Alternative splicing plays two well-established roles in gene 
regulation. First, it expands protein diversity by generating differ- 
ent mRNA isoforms from protein-coding genes. Second, alterna- 
tive splicing attenuates gene expression by producing substrates 
for nonsense-mediated decay (Ni et al. 2007; Pan et al. 2008; Wang 
et al. 2008). The work presented here suggests a third and un- 
expected consequence of alternative splicing — the translational 



control of resultant mRNA isoforms (AS-TC). Our data suggests 
that splicing generates transcripts with different affinities for the 
translational machinery. We hypothesize that this isoform-specific 
mRNA translation maybe regulated depending on the cellular needs. 

Translational control is predominantly achieved through 
regulation of initiation (Lackner and Bahler 2008). C/s-acting 
translational enhancers have been characterized in diverse species 
and influence either cap-dependent or independent translation 
(Pfeiffer et al. 2012). Eukaryotic viruses and many cellular mRNAs 
contain regulatory elements in 5' UTRs such as internal ribosome 
entry sites and tandem oligopyrimidine elements that regulate 
translation initiation under a variety of cellular conditions (Jeff eries 
et al. 1997; Jackson 2005). One possibility is that alternative 
splicing alters the topology of these c/s-acting regulatory elements, 
resulting in isoform-specific fates (Sanford et al. 2004; Staudacher 
et al. 2011). Another possibility is that multifunctional c/s-acting 
elements, such as the FN1 EDA splicing enhancer, may directly 
couple splicing decisions to translational control (Lu and Cullen 
2003; Nott et al. 2004; Sanford et al. 2004; Barberan-Soler et al. 
2009). There is also new evidence supporting the idea of functional 
microRNA targets sites within coding exons (Hausser et al. 2013). 
This model, if true, provides an elegant mechanism(s) for coupling 
alternative splicing with isoform-specific mRNA translation. An 
alternative model, not considered here, is that c/s-elements regu- 
lating isoform-specific translation function in auto-regulatory 
feedback loops as described for the TYMS and DHFR mRNAs (Lin 
et al. 2000; Tai et al. 2004). 

RNA secondary structures can also regulate translational 
control of mRNA isoforms (Bonnal et al. 2003; Goodarzi et al. 
2012). A prime example of this type of regulation occurs during the 
localization of the Oskar mRNA to the posterior pole of the Dro- 
sophila oocyte. Splicing of the first intron of the Oskar pre-mRNA 
generates a secondary structure that is required for precise locali- 
zation of the mRNA. This stem-loop structure encompasses the 
exon-exon junction and creates binding sites for a variety of trans- 
acting factors (Hachet and Ephrussi 2004; Yano et al. 2004; Ghosh 
et al. 2012). It is not known if similar coupled-regulatory mecha- 
nisms exist in mammalian genomes. It will be interesting to de- 
termine if local structural elements in alternatively spliced exons 
(both the "skipped" and "included" products) influence trans- 
lation initiation. 

We think our data support a "coupled fate hypothesis," sug- 
gesting that alternative splicing and/or mRNP assembly pathways 
in the nucleus can result in different fates for messages (Fig. 6); this 
may result in active translation, sequestration or mRNA decay 
(Dreyfuss et al. 2002). We propose that alternative splicing is not 
only important for expanding the protein coding potential of the 
human genome but also for altering the c/s-regulatory composi- 
tion of mRNAs leading to an unexpected level of post-transcrip- 
tional regulation in the cytoplasm. 

Methods 

Cell culture, fractionation, and polyribosome profiling 

Human embryonic kidney (HEK) 293 T cells were cultured in 
complete medium (DMEM 10% FBS) at 37°C, 5% CO z . Prior to 
lysis, cells were treated with cycloheximide (100 |xg/mL), 10 min at 
37°C 5%C0 2 . Cells were washed three times with ice cold PBS and 
lysed in ice cold buffer A (0.5% NP40, 20 mM Tris HC1 pH 7.5, 100 
mM KC1 and 10 mM MgCl 2 ). Lysates were passed three times 
through a 23 G needle and incubated on ice for 7 min. Extracts 
were then centrifuged at 10K rpm for 7 min at 4°C. The superna- 
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Figure 6. Alternative splicing generates mRNA isoforms with diverse cytoplasmic fates. Alternative 
splicing expands the protein coding capacity of human genes (upper panel). Alternative splicing can 
also be coupled to nonsense mediated decay (AS-NMD, middle panel). This mechanism can attenuate 
gene expression post-transcriptionally by producing unstable PTC-containing isoforms. In this paper, 
we present a novel mode of splicing-based gene regulation linking alternative splicing and translational 
control (AS-TC, lower panel) of the resultant mRNA isoforms. Our data suggest the intriguing hypothesis 
that alternative splicing of both coding and noncoding exons alters the c/'s- regulatory landscape of 
mRNA isoforms leading to different levels of polyribosome association. 



tant was collected as crude cytosolic extract. Cytosolic extracts 
were overlaid on 10%-50% sucrose gradients prepared in 20 mM 
Tris HC1 pH 7.5, 100 mM KC1 and 10 mM MgCl 2 buffer (prepared 
using the Gradient Station, Biocomp Instruments). Gradients were 
then ultracentrifuged at 40K rpm for lh 20 min at 4°C using an 
SW41 in a Beckman ultracentrifuge. Individual polyribosome 
fractions were subsequently purified using a Gradient Station 
(Biocomp Instruments). RNA was purified from both cytosolic 
and polyribosomal fractions using TRI Reagent LS (Sigma), 
converted to double stranded libraries using the TruSeq polyA + kit 
(Illumina) and sequenced using the Illumina HiSeq2000 platform. 

High-throughput data analysis 

High-throughput RNA sequencing was performed at the UCSC 
Genome Sequencing Center, producing —50-100M reads per 
sample (Supplemental Table 3). Mapping of single-end ~35-bp 
reads was performed using Bowtie and TopHat (Langmead et al. 
2009; Trapnell et al. 2009) to the hgl9 build of the human genome 
(February 2009 GRCh37, NCBI Build 37.1). Mapping parameters 
are listed in Supplemental Table 1. All data collection and parsing 
was performed with the Perl programming language and statistical 
analyses with the R statistical language, version 2.14.1. All external 
library packages used are available on CPAN or CRAN, respectively. 

Estimation of alternative splicing 

We used the Mixture of Isoforms (MISO) Bayesian Inference model 
vl.0 for quantification of alternative RNA processing events (Katz 



et al. 2010). The MISO algorithm counts 
the numbers of reads that are common 
to both events and the reads that are ex- 
clusive to one isoform or the other, in 
order to estimate the percent-spliced-in 
(^P) in a given sample. For a given RNA 
processing event, we calculated by 
subtracting ^ from two different samples. 
We employed an "event-centric" ap- 
proach, considering annotated alterna- 
tive mRNA processing events with at least 
20 exclusive reads and a maximum con- 
fidence interval width of 0.45 for the 
MISO estimated posterior distribution of 
percent-spliced-in ("Psi", Supplemen- 
tal Table 5 lists all of the MISO outputs). 
The MISO events database was down- 
loaded from the MISO website (http:// 
genes.mit.edu/burgelab/miso). Orienta- 
tion of ^ vs. 1-^ represents the included 
vs. the excluded isoform for alternative 
5' and 3' splice sites (A5SS, A3SS), skipped 
exons (SE), and retained introns (RI). 
For alternative first and last exons (AFE, 
ALE), ^ and 1-^ refer to the proximal 
and distal exon sets, respectively. In the 
case of the analysis of functional ele- 
ments in Figure 5, the orientation of 
AFE, ALE, and MXE events may be 
switched such that the functional ele- 
ment in question is contained within 
the ^ isoform rather than 1-"*P. Classifi- 
cation into 5 ' UTR, CDS, and 3 ' UTR events 
was based upon the alternative region's 
orientation with respect to the March 
18th, 2012, update of human NCBI Con- 
sensus CDS database (NCBI Build 37.3) (ftp://ftp.ncbi.nlm.nih. 
gov/pub/CCDS/). 

False discovery rate calculations 

To empirically estimate the FDR based on biological repro- 
ducibility, we considered the filtered events common to both 
replicates. For each value of a: in 0.01 to 0.3 increasing by 0.01, we 
perform 100 bootstrapped label permutations of replicate 2. For 
each of the 100 shuffled sets, we determined the number of 
events where both replicates have at least a BF > 5, and|A^| > x, and 
following the same direction of change. The FDR for each of the 
100 permutations of each value of x was then estimated as FP/(FP + 
TP), which is the number of permuted events passing the criterion 
divided by the number of real events passing the criterion. This 
was plotted against the number of true positives, the number of 
real events passing the criterion minus the number of permuted 
events. 



Reverse transcription-polymerase chain reaction 

RNA was purified from cytoplasmic extracts and polyribosome 
fractions using TRI Reagent LS (Sigma) following manufacturer's 
instructions. cDNA synthesis was performed on each fraction 
using random hexamers (Applied Biosystems) following manu- 
facturer's instructions. PCR was performed using Taq Titanium 
(Clontech). For cytoplasm to polyribosme comparisons, rate of 
inclusion was calculated using the following formula: %Inclusion= 
Included/ (Skipped+Included) in nmol as determined using an 
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Agilent 2001 Bioanalyzer. Statistics were done using a Mann- 
Whitney one-tailed t-test. 

Oligos 

MOV10 FOR: GGCGACTTTCAGTTTCATTTCCAC 
MOV10 REV: TCTTGGCTAGCAGTGACTTGTGG 
THOC5 FOR: GCTCCGGGACTGAGGTGCT 
THOC5 REV: CTCGGTGTCAGATCGATTCCG 
DHX9 FOR: CAGAAGAAGTGGATTTAAATGCTGG 
DHX9 REV: CCCAGCTGCTTGATATAAATGGTC 
KARS FOR: GGAAGATGGCGGCCGTGCA 
KARS REV: CTCTTCCTCAGGACCCACAC 
HNRNPR FOR: TGCTGCCGGGATCGCGTCT 
HNRNPR REV: AGTCACTTTCCTTGAACTGCTGTAG 
SGOL1 FOR: GCAAATTCCTCTTGAAGAAACTGAAC 
5GOL2 REV: AGGGTGGGCTCCTTATAGTTCAC 
NDUFS1 FOR: ATATGTTAAGGATACCTGTAAGAAAGGC 
NDUFS1 REV: GGGCTTTCTCAATTTCAACAAGGC 
DDX50 FOR: AATGCCTGGGAAACTCCTCT 
DDX50 REV: TGGTGCCTTGACTTCCTTCT 
DGUOK FOR: CTCTCCATCGAAGGCAACAT 
DGUOK REV: TTTCCAAGACTTTGGGCAGT 



Data access 

The Frac-seq data from this study have been submitted to the NCBI 
Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/ 
geo/) under accession number GSE45237. 
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